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^^ ' I discuss the sensitivity of future medium baseline reactor antineutrino experiments on the neu- 

'_ ' I trino mass hierarchy. By using the standard x^ analysis, we find that the sensitivity depends strongly 

'"— ' . on the baseline length L and the energy resolution {5E/E)^ = (a/^E/MeV)'^ + 6^, where a and b 

parameterize the statistical and systematic uncertainties, respectively. The optimal length is found 



(N 



Jh ' to be L ~ 40 — 55 km, the larger resolution the shorter optimal L. For a 5 kton detector (with 12% 



Ph' 



o 



^ 



weight fraction of free proton) placed at L ~ 50 km away from a 20 GWth reactor, an experiment 



■^^ ' would determine the mass hierarchy with (Ax^)min ~ 9 on average after 5 (15) or more years of 



running with the (a, 6) = (2,0.5)% ((3,0.5)%) energy resolution. This type of experiment can also 
measure the relevant mixing parameters with the accuracy of ~ 0.5%. 



^ ! I. INTRODUCTION 

Ph. 

Ph' Now that a large Q\z has been measured at Daya Bay [l|, Q and RENO [3[ experiments accurately, neutrino 

2;^ \ physics enters a new era. One of the next challenges is determination of the mass hierarchy. 

Among many ideas proposed, the medium baseline reactor antineutrino experiment [4l4lCll| has stimulated 
various re-evaluations of its physics potential and sensitivity recently. Some works utilize the Fourier transform 

CN , technique jll| - [l3| . first discussed in refs. 0-13 j to distinguish the mass hierarchy. The main advantage of this 

^ ■ technique is that the mass hierarchy can be determined without precise knowledge of the reactor antineutrino 

^X \ spectrum, the absolute value of the large mass-squared difference |Am§;^|, and the energy scale of a detector. 

I^ ■ Although interesting and attractive, this technique is somewhat subtle to incorporate the uncertainties of the 

ly— , \ mixing parameters and to estimate its sensitivity to the mass hierarchy. 

• . On the other hand, some works adopt the x^ analysis [l3, [l3l - [l5l | and new measure based on Baycsian 

iSs ' approach [16|. These methods utilize all available information from experiments, and it is straightforward to 

-Y-\ . incorporate the uncertainties to evaluate the sensitivity, providing robust and complementary results to the 

,_^ ' Fourier technique. 

• • \ In this proceedings, we analyze the sensitivity of medium baseline reactor antineutrino experiments to the 

T'"! , mass hierarchy for the baseline length of 10-100 km and the energy resolution i^EjE)'^ = la/^/E/McVj +6^ 

5_] ] in the range 2% < a < 6% and b < 1% with the x^ analysis. The optimal baseline length and the expected 
Cd ■ statistical uncertainties of the neutrino parameters, sin^ 20i2,sin^ 26'i3, ^"^ii ^^-^ Arn^j^, are also estimated. 

II. REACTOR ANTINEUTRINO FLUX 

In this section, we briefly discuss the evaluation of how many electron antincutrinos, Pg, would be detected 
at a far detector with a medium baseline length from a reactor. 

In a nuclear reactor, antincutrinos are mainly produced via beta decay of the fission products of the four 
radio-active isotopes, ^^^t/, ^^^t/, '^^^Pu and ^^^P„, in the fuel. The flux of antincutrinos with energy E,^ (MeV) 
at a reactor of P (GWth) thermal power is then expressed as Q 

^^^^^HE.)x6.2ixlO^\ (1) 

where /;; and e^ are the relative flssion contribution and the released en ergy per flssion of the isotope fc. 



respectively [17[. 4){E^) is the number of antincutrinos produced per flssion |18 | 
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This rate is then modulated by neutrino osciUation. The De survival probability is expressed as 

3 2 



^^ Uei exp ( -i^E' I ^*' 

= 1 - cos^ 6*13 sin^ 26'i2 sin^ (A21) 
-sin2 20i3sin2(|A3i|) 
- sin^ 6*12 sin^ 26*13 sin^ (A21) cos (2| A31I) 

± !!!!^ sin^ 2013 sin (2A21) sin (2| A31 1) , (2) 

where Uei is the neutrino mixing- matrix clement relating the electron neutrino to the mass eigenstate Vi. The 
variables mi and Ei arc the mass and energy of the corresponding mass eigenstate, while 9ij represent the 
neutrino mixing angles. The oscillation phases Aij are defined as 



Am|L 
4^„ 



^'^^^-n^^ (Arr4^m?-m^) (3) 



with a baseline length L. The plus or minus sign in the last term of eq. ([5]) corresponds to the mass hierarchy; 
the plus sign is for normal hierarchy (NH), and the minus sign for inverted hierarchy (IH). Note that this last 
term is the only source of the mass hierarchy difference. We have neglected the matter effect because it is small 
enough for the energy range and the baseline lengths we concern in this study [l9|. 

Similar as the current reactor experiments, such as Daya Bay [l|, y], RENO Q and Double Chooz [20], 
future medium-baseline reactor-antineutrino experiments can also use free protons as targets to detect electron 
antineutrinos via the inverse neutron-beta-decay (IBD) process, producing a neutron and a positron. The 
threshold neutrino energy of this process is 

Ethi ^ rn„ - nip + me. (4) 

The produced positron then interacts with scintillator, converting its kinetic energy to photons. Eventually, 
the positron annihilates with an electron in the detector and emits two 0.5 McV photons. The energies of 
photons are accumulated as the visible energy, -Evis, which is the sum of the positron's total energy and an 
electron's rest energy, 

E^i^ ^Ee -f me - {E^ - 0.8) MeV. (5) 

The observed antincutrino distribution by a detector with Np free protons after an exposure time T can then 
be expressed as 

dN N T f°° dN 

dE°P^ ^^ J '^^''1^/''^^' E,)<jisd{E,.) G(E, - O.SMcV - K?/, SE^,,), (6) 

where aiBDiE^) is the cross section of the IBD process [2l|, G is the detector response function with the energy 
resolution SE^i^, and E°f^ is the observed visible energy by the detector. In this study, we take the normalized 
gaussian function as the response function, i.e., 

CiE„. - E-..E.,.) . -^e.p{-(^||^ j ^ (7) 

The detector energy resolution, SE^is, is parameterized into two parts. 



E^ 



\ \yE~/McV^ 



b'- (8) 



The first term in the square root represents the statistical uncertainty, and the second one gives the systematic 
uncertainty |22| . 
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TABLE I: The input values yi"P"' and their uncertainties SY taken from refs. [3, 0, [13]. The uncertainty of sin^ 2013 
can be 5% or less after 3 years running of Daya Bay experiment [2J| . 



Y 


sin^ 26ii2 


sin^ 26li3 


Amii eV" 


lAmiileV^ fsys 


\/-input 

5Y 


0.857 
0.024 


0.089 
0.005 


7.50 X lO^"^ 
0.20 X lO"-"^ 


2.32 X 10"^ 1 
0.1 X 10"^ 0.03 



III. THE SENSITIVITY TO THE MASS HIERARCHY 

After obtaining the energy distribution of reactor antineutrinos^_we would like to estimate the sensitivity of 
determining the mass hierarchy using the standard x^ analysis [H, H, uQ, [l3l - [l5| . 
To set the stage, we introduce the x^ function as 

A Apara ' Astat" \^/ 

The first term summarizes the prior knowledge on fitting parameters. In reactor antineutrino experiments, 
these are the mixing angles, sin 29i2 and sin 20i3, and the two mass-square differences, Arn^i and |Ato|]^|. In 
this study we also consider the event-number normalization factor /sys, assuming the 3% uncertainty. Their 
contributions look like, 

2 r (sin^ 26*12) ^^ - (sin^ 26'i2) '"?"* 1 ^ f (sin^ 26*13) ^* - (sin^ 26*13 



^P"'" 1 Ssm^26i2 J I S sin' 291 



3 



(Amli)fi'-(AmiJ-P"'\' ^ f {\Aml,\)^' ~ i\Aml,\y-P-'\^ , ( fj^~fj^ 



(S|Am|J J ' \ Sfsys 



(10) 



The input values y™?^* and their uncertainties 6Y are listed in Table HI 

The second term in ^ represents the statistical fluctuation. When we introduce binning w.r.t. E°f^, it looks 
hke 



jV.fi* - 7V^H(IH) 



with the summation running over all the bins. Here, N^ is the event number for the zth bin when the 

hierarchy is NH (IH), while A^f* is the theoretical prediction of the event number cither with right or wrong 
mass hierarchy, calculated as a function of the four model parameters and the normalization factor /sys, which 

are all varied under the constraints of ([TO)). In this study we prepare the data A^^ ^ -^ by using eq. ([6]) with 
the input values of the five parameters for each mass hierarchy. 

In the limit of infinitely many events, the bin size can be reduced to zero, and the sum ()11|) can be replaced 
by an integral. Although a finite bin size is required for actual experiments, we adopt this zcro-bin-sizc limit as 
a measure of the maximum sensitivity. 

We then define A^^ as 

Ax' = X' - Xmin, (12) 

where x^;„ is the minimum of x^ 1 which is obviously zero in our approximation of neglecting statistical fiuctu- 
ations in data, N,- . When wrong mass hierarchy is assumed in the fit, the minimum of Ax^, (Ax^)ijiin, 



will deviate from zero, and the wrong mass hierarchy can be rejected with the significance ^y(Ax 



2\ . 
mm- 
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IV. RESULTS 

In this section, we discuss the sensitivity to the mass hierarchy, the optimal length and the statistical uncer- 
tainties of the neutrino parameters. All our results are obtained by assuming a reactor of 20 GWth thermal 
power, a far detector of 5 kton fiducial volume with 12% weight fraction of free proton and 5 years exposure 
time. 




E^ [MeV] 



E„ [MeV] 



FIG. 1: The energy distribution of reactor antineutrinos with 20GWth-5kt (12% free-proton weight fraction) -Syrs ex- 
posure and the baseline length L = 30 km (left) and 50 km (right). Upper: The case with exact E^, measurement, 
where the dashed blue and dashed red curves are for NH and IH, respectively. The solid curve shows the best fit of IH 
assumption to the NH data. The red arrow points out the energy at which the difference due to mass hierarchy vanishes. 
Lower: 6/V-Bvis % energy resolution case. 



Figures [T] shows energy distributions for L = 30 km (left) and 50 km (right), in which the exact E^ measure- 
ment is assumed in the upper panel, whereas in the lower panel the energy resolution of a = 6% with 6 = 
in eq. (|8]) is assumed. The dashed blue curve corresponds to the NH case, and the dashed red curve to the IH 
case, while the solid curve is obtained using the parameter values fitted to the NH data with the "wrong" IH 
assumption. The red arrow points out the energy at which the difference due to mass hierarchy vanishes (the 
degeneracy point). 

At L = 30 km, the solid curve almost coincides with the dashed blue one even with the exact energy 
measurement, implying that it is almost impossible to distinguish the mass hierarchy by experiments at L = 30 
km. This is because the small phase shift between the NH and IH predictions can be absorbed by a small shift 
in |Am§i| by a fraction of its present uncertainty, 0.1 x 10^'^eV^. 

The situation changes when the second peak of the mass-hierarchy dependent term appears in the energy 
range. The mass hierarchy difference can no longer be absorbed by a shift in |A?n|j^| since the relative phase 
difference between the NH and IH oscillations changes across the degeneracy point. There is no way to make 
the differences on the both sides compensated, resulting in the distinct mismatch between the dashed blue curve 
(the NH data) and the solid curve (the best fit under the IH assumption) as shown in the upper panel of the 
right plot in Fig. [1] where the antineutrino energy is exactly measured. 

Once the finite energy resolution is introduced, the phase difference in the lower energy side of the degeneracy 
point is significantly smeared out as it oscillates faster w.r.t. E^, at the low energy. Hence it is easier for one 
oscillation period to be covered by a sizable Gaussian profile of the detector response function. The remaining 
difference in the higher energy side can then be absorbed by a small shift in jAjTigj^l, resulting in an excellent 
fit (solid curve) to the NH data (blue dashed curve) in the lower panel of the right plot in Fig. [TJ shown for 
6%/\/E/MeV energy resolution. 

The left plot in Fig. [2] shows the resulted (Ax^)min value as a function of the baseline length L, for several 
energy resolutions, a = 2,3,4, 5 and 6% (with 6 = 0) in eq. ([5|), from the top to the bottom. Solid curves are 
for NH, while dashed curves are for IH. The results clearly show that t he mass hierarchy can be determined by 
those experiments only if the energy resolution of the detector is 3%/y^E/McV or better, and that the optimal 
baseline length (as shown by the cross symbol) is around 50 km for that resolution. The small (Ax^)min for the 
baseline length L < 40 km and L > 80 km is due to a shift in |A77i§]^| and low statistics, respectively. For the 
a = 5 and 6% cases (Ax^)niin stays almost zero at all L. 
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FIG. 2: (Ax^)min for mass hierarchy discrimination shown as a function of the baseline length L for 20GWth-5kt (12% 
free-proton weight fraction) -Syrs exposure. The solid curves are for the NH cases and dashed curves for the IH cases. The 
cross symbols mark the optimal baseline lengths. Left: The energy resolution in eq. ([8} is varied with a = 2 to 6% and 
6 = 0, from the top to the bottom. Right: The energy resolution is varied with a = 2% and h — 0%, 0, 5%, 0.75%, 1%, 
from the top to the bottom. 



The right plot in Fig. [5] shows the (Ax^)min value as a function of the baseline length L for different b 
values with a = 2%. The curves from the top to the bottom are obtained for b = 0%, 0.5%, 0.75% and 1%, 
respectively. The effect of the systematic uncertainty is significant as discussed in ref. |13| , reducing the peak 
value of (Ax2)mm from 11.0 {b = 0) to 9.7 {b = 0.5%), 8.4 {b = 0.75%) and 6.9 {b = 1%) for NH. The optimal 
L shortens from 51 km for (a, b) = (2, 0)% to 47 km for (a, b) = (2, 1)%. 

In addition, the neutrino parameters, sin 2^12, Attiji and | A77i§]^ |, can be measured accurately with statistical 
uncertainties shown in Fig. [S] We find 
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FIG. 3: The statistical uncertainties of the neutrino model parameters measured by this experiment as functions of the 
baseline length L after 20GWth-5kt (12% free-proton weight fraction) -Syrs exposure. The results are shown for both 
hierarchy (NH by solid and IH by dashed curves) and for the energy resolution of eq. ((8| with (a, 6) = (3, 0.5), (3, 1) and 
(6, 1)%. 



(5sin2 26li2 -4 x 10"^ (0.5%), (13a) 

,5Am^i -3 X lO^^cV^ (0.4%), (13b) 

<5|Ato^i| -7 X lO^^cV^ (0.3%), (13c) 

with the energy resolution of (a, 6) = (3, 0.5)% at L = 50 km; the percentage values in the parentheses denote 
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the relative accuracy of the measurement. Those uncertainties are almost independent of the mass hierarchy 
and of the energy resolution, with the only exception of the |A77i|]^| uncertainty for which the larger resolution 
results in the larger uncertainty, sin 20i3 and |A77i§]^| are measured most accurately around L ^ 1 km. 

V. DISCUSSIONS AND CONCLUSION 

In this proceedings we have investigated the sensitivity of medium-baseline reactor-electron-antineutrino 
oscillation experiments for determining the neutrino mass hierarchy by performing the standard x^ analysis. 

We have carefully studied the impacts of the energy resolution {SE/E)^ ~ ( a/^i?/McV] + b^ and find that 

the sensitivity strongly depends on it. The optimal baseline length is found to depend slightly on the energy 
resolution, preferring the length slightly shorter than 50 km for the energy resolution of (a, 6) = (2, 0.75)% and 
(2, 1)%. At the optimal baseline length, the energy resolution better than the 3%/y'i?/MeV level is needed to 
determine the neutrino mass hierarchy pattern. For a 5 kton detector (with 12% weight fraction of free proton) 
placed at L ~ 50 km away from a 20 GWth reactor, an experiment would determine the mass hierarchy with 
Ax^ ~ 9 on average after five or more years of running if the energy resolution of (a, 6) = (2, 0.5)% is achieved, 
while a factor of three larger or longer experiment is needed to achieve the same goal for the energy resolution 
of (a,&) = (3,0.5)%. 

It is also found that this experiment can measure the neutrino parameters, sin 20i2, ^rn^i and IAttij]^!, very 
accurately as shown in (fT3|) for an experiment of 20 GWth-5kt (12% free-proton weight fraction) •5yrs at L ~ 50 
km. 
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